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Abstract 

We determine the crossover exponents associated with the traceless tenso- 
rial quadratic field, the third- and fourth-harmonic operators for O(n) vector 
models by re-analyzing the existing six-loop fixed dimension series with the 
pseudo-e expansion. With this approach we obtain accurate theoretical esti- 
mates that are in optimum agreement with other theoretical and experimental 
results. 
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I. INTRODUCTION 



For many years the critical behavior of 0(n) vector models has caught a lot of attentions 
since most of the physical systems undergoing second-order phase transitions belong to the 
O(n) universality classes (see Ref. [1] for a recent review). Thus a precise determination of 
universal quantities as critical exponents and amplitude ratios has became necessary. The 
critical behavior of these physical systems may be obtained by field theoretical investigations 
based on the Landau-Ginzburg- Wilson Hamiltonian 



H = Jd d x ^dj- dj+ l -r$- + • 0) 2 



(1) 



where <f>(x) is a n-component real field. An interesting issue is to determine how the critical 
properties are influenced by the addition of a perturbation term to the Hamiltonian (1) 



n 



J d d x [^0 • + 0) 2 + h p v\ , (2) 



where h p (x) is an external field coupled to V(x). In fact, if V is an eigenoperator of the 
Renormalization Group (RG) transformations, the singular part of the Gibbs free energy 
becomes a scaling function in the limit of reduced temperature t — > and h p — > 0, and can 
be written as 

F sing (t,h p ) w \t\ d ^ (h p \t\-^) , (3) 

where P = y p v is the crossover exponent associated with the perturbation V and y p is the 
RG dimension of V. Moreover, one usually defines the indices f3 p and 7 P which describe 
the low-temperature singular behavior of the average (V(x)) ~ \t\P" and of the suscettivity 
Xv — I d d x(V(x)V(0)) c ~ t~~ Jp . They satisfy the scaling relations 

P P = 2 - a - <f> p , 7 P = -2 + a + 20 p . (4) 

Among the perturbation operators, particularly important from the experimental and 
phenomenological point of view are the so-called harmonic ones [2-4] 

V 2 (x) = (f)i(x)(t)j{x) - Sij - $(x) ■ ${x), 

71/ 

V 3 (x) = (fra<frb<frc ~ ^—^(<Pa5bc + (frb^ac + <PJab), (5) 

Vi{x) = 4>a4>b(t>c<t>d ~ ^—^7{Sab<Pc<Pd + 5perm.) + - — ^rrt- : oA {S a bS c d + 2perm.) , 
n + 4 (n + 4)(n + 2) 

called second, third, and fourth harmonic operator respectively. In the following we will 
denote the crossover exponent of Vi as 0j and its RG dimension as yi. Higher order harmonic 
operators are generally reputed to be irrelevant at the three-dimensional O(n) fixed point 
[3], thus we will not consider them here. 

The crossover exponent 2 associated with the traceless tensor field Vi{x) reveals the in- 
stability of the 0(n)-symmetric theory against anisotropy [2,5-7]. It characterizes the phase 
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diagram at the multicritical point where two critical lines O(n) and O(m) symmetric meet. 
In some cases this gives rise to a critical theory with enlarged 0(n + m) symmetry [8-11]. 
Multicritical behavior arises in several different contexts in physics: in anisotropic antiferro- 
magnets in a uniform magnetic field [6,11], in high T c superconductors (see e.g. Ref. [12] and 
note that in the SO (5) theory of superconductivity [13] the multicritical point is effectively 
0(5) symmetric), in colossal magnetoresistance materials [14], in certain theories of strong 
interactions [15], etc. This list is far from being exhaustive, we only quote some examples. 
For the XY model (n = 2) the traceless tensor field V 2 (x) and its correlation function 
are connected with the second-harmonic order parameter in density- wave systems [16,17], 
which characterizes some liquid crystals at the nematic-smectic-A transition [16-23]. The 
structure factor of the secondary order parameter V 2 , that within RG methods has been 
determined in- [16,24] and out-of-equilibrium [25], has been experimentally measured using 
X-ray scattering techniques [22,23]. Finally the RG dimension y 2 enters in the study of 
crossover effects in diluted Ising antiferromagnets with n-fold degenerate ground state [26], 
in models with random anisotropy [27], at certain quantum phase transitions [28], and in 
other more complicated situations [29]. Even this list is far from being exhaustive. 

The third-harmonic crossover exponent determines the phase diagrams at the smectic- 
A hexatic-B point in liquid crystals [18,23], in materials exhibiting structural normal- 
incommensurate phase transitions [30-32], and at the trigonal-to-pseudotetragonal tran- 
sition [33]. For n = 0, 03 is related to the partition function exponent of nonuniform star 
polymers with three arms [34]. Finally, for n > 2 it determines the stability of O(n) fixed 
points against n + 1-state Potts-like perturbations [35] as it happens in the presence of stress 
or particular magnetic fields [33,36]. 

The fourth-harmonic exponent 04 is mainly related to the stability of the O(n) fixed 
point against fourth-order anisotropy [8], as e.g. the cubic one [7]. It is worth mentioning 
that for n = 1, even if the operators Vi(x) can not be defined through Eqs. (5), all the 0j 
have non trivial values. This fact has an interpretation in terms of a gas of n-colors loops 
(see e.g. [37]) in the limit n — > 1. 

The exponents 0j and yi with i — 2,3,4 have been analyzed in the past with different 
theoretical methods, in the framework of the e-expansion [38-41,2,3,18,8,34], from the anal- 
ysis of high-temperature expansion [9] , by means of Monte Carlo simulations [42-44] , in the 
large n approach [3,45,46], and in the fixed- dimension perturbative expansion [18,24,34]. 

The aim of this paper is to determine the crossover exponents 0j and the RG dimensions yi 
by re-analyzing the three-dimensional six-loop perturbative series [24,34] with the pseudo- 
e expansion trick [47], since in many cases this method provided accurate results in the 
determination of critical quantities (see, e.g., Refs. [48-52]). The idea behind this trick is 
very simple: one has to multiply the linear term of the j3 function by a parameter r, find 
the fixed points (i.e. the zeros of the (3 function) as series in r and analyze the results as in 
the e expansion. The critical exponents are obtained as series in r inserting the fixed-point 
expansion in the appropriate RG functions. With this trick the cumulation of the errors 
coming from the non-exact knowledge of the fixed point and from the uncertainty in the 
resummation of the exponents is avoided. The obtained pseudo-e series are believed to be 
asymptotic, so an appropriate resummation is usually needed to have reliable estimates. [48] 
However, it often happens that up to six-loop order, the pseudo-e expansions do not yet show 
their asymptotic nature, rendering effective a non-resummed evaluation based on simple 
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Pade. 

The paper is organized as follows. In Sec. II we analyze the quadratic crossover expo- 
nents, in Sec. Ill the cubic and quartic ones. In Sec. IV we report all the pseudo-e estimates 
for harmonic exponents and compare them with other theoretical and experimental ones. 



II. QUADRATIC CROSSOVER EXPONENTS 

The six-loop RG perturbative series in the three-dimensional approach for the second 
harmonic operators were computed in Ref. [24], whereas the /3 function (necessary to find 
the stable fixed point) is reported for general n in Ref. [53]. By using these series, one 
obtains the pseudo-e expansion of all the second-harmonic exponents. 

We first consider the crossover exponent 02- As a typical example, the perturbative 
expression in the parameter r for n = 2 reads 

r 317t 2 

0™= 2 = l + _ + + 0.01688r 3 + 0.01005r 4 + 0.00227r 5 + 0.00457r 6 + 0(r 7 ) . (6) 

2 10 6750 

At least up to the presented number of loops the series does not behave as asymptotic with 
factorial growth of coefficients. Although the series has not alternating signs, which is a key 
point to ensure some kind of convergence, one can try to apply a simple Pade summation. 
The results for n = 2 are displayed in Table I. All the approximants possess poles on the 
real positive axis. Some of them are close to r = 1 and the estimate of 2 on their basis 
should be considered unreliable. Anyway some of these approximants have poles "far" from 
t = 1, where the series must be evaluated. Thus one may expect the presence of such poles 
not to influence the approximant at r = 1. Indeed all such Pade results are very close since 
lower orders. Hereafter we choose as final estimate the average of those six-loop order Pade 
without poles in < r < 2, and as error bar we take the maximum deviation of the final 
estimate from the four-, five- and six-loop Pade. The five-loop estimates are analogously 
obtained, considering the maximum deviations up to three loops. Within this procedure we 
obtain 02 =2 = 1-178(15) at five-loop and 1.181(7) at six-loop. Although in good agreement 



TABLE I. Pade table for <^= 2 m pseudo-e expansion. The two integer numbers N and M 
denote the corresponding [N/M] Pade. The location of the positive real pole closest to the origin 
is reported in brackets. 





N = 


N = 1 


N = 2 


N = 3 


N = 4 


N = 5 N = 6 


M = 





1 


1.1 


1.14696 


1.16385 


1.17390 


1.17616 1.18073 


M = 


1 


1.11111 [10 ] 


1.18855(0!] 


1.1 7332^9 s] 


1.18868 [1 7] 


1.17683 M 


1.17166 [0 . 5] 


M = 


2 


1. 15870 M nj 


1. 17369 [o 7] 


1. 17925 [o 3] 


1.17952 [0 3] 


1.18113 [2 .i] 




M = 


3 


1.17021 [3 . 2] 


1. 19549 p 4] 


1.17952 [2 . 3] 


1.17920(0.!] 






M = 


4 


1.17818[o 7] 


1.17773[ot] 


1.18127 [2.03] 








M = 


5 


1.17771 [2 . 7] 


1.17814[o 7] 










M = 


6 


1. 18271 [o ,] 
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TABLE II. Pade table for y 2 =2 in pseudo-e expansion. The two integer numbers N and M 
denote the corresponding [N/M] Pade. The location of the positive real pole closest to the origin 
is reported in brackets. The final estimate is y~2 =2 = 1.763(4). 





N = 


N = 1 


N = 2 


N = 3 


TV = 4 


N = 5 


N = 6 


M = 





2 


1.8 


1.76089 


1.76705 


1.76215 


1.76707 


1.76041 


M = 


1 


1.81818 


1.75138(5.!] 


1.76621 


1.76432 


1.76461 


1.76424 




M = 


2 


1.77061 


1.76756 


1.76376[ 98 ] 


1.76459 


1.76443 






M = 


3 


1.76774 


1.77550 [0 . 6] 


1.76471 


1.76433 [13 ] 








M = 


4 


1.76322 


1.76505 


1.76405(3.4] 










M = 


5 


1.76630^ 


1.76444 












M = 


6 


1.76162 















with other theoretical estimates, we do not retain safe such estimates, since the presence of 
so many poles in the Pade can cause systematic deviations from the actual value. This can 
be traced back to the fact that the series has not alternating signs. To improve the estimates, 
one can try to resum the series by means of the Pade-Borel-Leroy (PBL) method (see e.g. 
[48]) or more advanced ones, but the monotonic character of the signs of the series makes 
the majority of the approximants to be defective. The resulting few good approximants do 
not allow a safe determination of the quantities analyzed. The same scenario is found for 
all other values of n. 

To achieve a reliable estimate of 2 , one has to consider series that have alternating signs. 
This can be done by considering the RG dimension y 2 = fyijv. The pseudo-e expansion of 
2/ 2 for general n is 

2 2(-392 - 78n + 5n 2 ) 2 

Vo = 2 r H — -r 

y2 (n + 8) 27(8 + nf 

362.271 + 65.4612 n + 16.8274 n 2 + 6.59007 n 3 + 0. 192062 n 4 , 
H t 

(8 + n) 5 

19417.6+10881.8 n+2401. 97 n 2 -29.3771 n 3 -83.6197 n 4 -5. 54754 n 5 -0. 0986411 n 6 4 

(n+8) 7 ' 

, 1.66461 10 6 +981069. n+215076. n 2 +30474.1 ra 3 +8065.54 n 4 +1734.33 n 5 +132.809 n 6 +4.22169 n 7 +0. 0573847 n g ^5 
^ (8+n) 9 ' 

" 1.54538 10 8 +1. 25687 10 8 n+4.50851 10 7 n 2 +9. 02640 10 6 n 3 +7.98957 10 5 n 4 -84675.5n 5 -31129.6n 6 -2974.43n 7 

(8+n) 11 

-133.244n 8 -3.26010n 9 -0.0369419n 
(8+n) 11 

that has alternating signs for n < 6. In fact, we get a small number of Pade with poles on 
the real positive axis, as one may appreciate from Table II where the results for n = 2 are 
displayed. The goodness of the Pade persists increasing n up to n ~ 6 while for higher n the 
results get worse. All the final data are shown in Table IV, where we also report the result 
for n = 16, that has to be taken with care since in this case the series is not alternating in 
signs. 

We resum the perturbative series by the PBL method too. The number of defective 
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approximants is very low and one may obtain a different estimate of the quantity y 2 ■ We 
consider the four-, five-, and six-loop approximants. The final estimate is the mean value 
between the maximum (M) and the minimum (m) value found. The error is (M — m)/2. 
During this average procedure, we discard the approximants that seem to be pathological. 
The results found with this method (displayed in Table IV) are very stable and always 
compatible with the Pade ones, but they have smaller errors. Only for n = 16 we are not 
able to give a PBL final result since most of the approximants are defective. This great 
stability within the PBL resummation makes us decide to report as final estimates (Table 
V) the simple Pade ones, in order to avoid underestimation of the uncertainties. 

Exploiting the scaling relations (4), we can apply the previous procedure to characterize 
the critical exponents f3 2 and 72. Unfortunately these series have no alternating signs for all 
values of n, resulting in a bad determination of their actual value. In Table IV we display 
2 , fl 2 , and 72 as obtained by using scaling relations from y 2 and the most accurate estimates 
of standard critical exponents in the O(n) universality class [54]. 



III. THIRD- AND FOURTH-HARMONIC CROSSOVER EXPONENTS 

In this section we consider the critical exponents of the third- and fourth-harmonic 
operators. 

The six- loop three dimensional series relevant for P 3 were calculated in Ref. [34]. Even 
in this case only the direct estimate of y 3 gives rise to a reliable result, since the series for 
/3 3 , 73, and 3 have not alternating signs. The pseudo-e expansion of y 3 for general n is 

3 6 144 + 44n - 14n 2 2 

Vs ~ 2 ~ ^+8 T 9(n + 8) 3 T 

404.981 + 281.073 n + 78.7256 n 2 + 20.9419 n 3 + 0.613212 n 4 - 
H T 

(n + 8) 5 

1 -22461.1-17805.3 n-3307.73 ra 2 +417.392 n 3 +244.052 n 4 +15.6783 n 5 +0.289476 n 6 4 
~T (n+8) 7 T 

I 1.67818 10 6 +1.44153 10 6 n+420045. n 2 +92283.5 n 3 +24863.3 n 4 +4900.83 n 5 +366.716 n (i +11.7112 n 7 +0. 162327 re 8 5 
"I (rt+8) 9 ' 



TABLE III. Pade table for y 3 ~° in pseudo-e expansion. The location of the positive real pole 
closest to the origin is reported in brackets. 
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N = 1 


N = 2 


N = 3 


N = 4 


N = 5 


N = 6 


M = 





1.5 


0.75 


0.71875 


0.73111 


0.72040 


0.73290 


0.71383 


M = 


1 


1 


0.71739 [24 ] 


0.72761 


0.72537 


0.72617 


0.72535 




M = 


2 


0.84706 


0.73152 


0.72529(47] 


0.72603 


0.72573 






M = 


3 


0.78599 


0.71920 [3 . 5] 


0.72621 


0.72572 [60] 








M = 


4 


0.75533 


0.73227 


0.72530 [6 . 5] 










M = 


5 


0.74240 


0.68973[i. 3] 










M = 


6 


0.73215 
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Tn 0.016 
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a, 

>* 0.012 
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[3/1] 
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- [1/3] 
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[1/5] 
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1 1 1 ' 1 

12: 



FIG. 1. (Color on-line) Non-defective four-, five, and six-loop PBL approximants for the RG 
dimension 2/4 for n = 2,3. 



1.63843 10 8 +1. 78385 10 8 n+7.35887 10 7 n 2 +l. 53121 10 7 n 3 +968499.n 4 -318541.n 5 -87011.6n 6 -8034.60ra 7 

(n+8) 11 

-360.085n 8 -8.93464n 9 -0.102502n 10 " 



(n+8) 1 



(8) 



which has alternating signs for n < 5. To show the goodness of the Pade summation we 
report in Table III the data for n = 0. The final estimate from this table is y% = 0.725(29) 
[y 3 = 0.731(35)] at six-loop [five-loop]. These estimates are obtained using the procedure 
outlined in the previous section. Similar good Pade tables are found for higher values of n, 
up to n ~ 5. All the final results are reported in Table IV, where we also show the PBL 
ones for a comparison. Again the PBL uncertainty seems too small to be considered safe. 

It is worth noting that for the partition function exponent p of non-uniform star polymers 
with three arms [34], we obtain the pseudo-e series 

p = 3(7 + u)/2 + 03 = 3 + 0.009278 r 3 + 0.005889 r 4 - 0.014350 r 5 + 0.015458 r 6 , (9) 

16 

which (by means of simple Pade) leads to p = 3.055(11). This value compares well with 
other estimates [34] and with the one obtained from 2/3 and the most accurate theoretical 
estimates of 7 and v [54] leading to p = 3.043(18). 

Finally, let us consider the fourth-harmonic exponent. It has been shown in Refs. [55,8] 
that the RG dimension t/ 4 is related to the exponent characterizing the stability of the 0(n) 
fixed point against cubic anisotropy. Thus we can use the six-loop series of the cubic model 
[56] to obtain the pseudo-e expansion 

12 4(680 + 62n + 23n 2 ) 2 
y* ~ 1 - ^ T + 27(8 + n)3 T 

-5200.56 - 777.127 n + 33.2649 n 2 + 35.7448 n 3 + 1.25111 n A , 
H t 

(n + 8) 5 
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TABLE IV. Final estimates from six-loop pseudo-e expansion. The values of j/j are calculated 
directly from the series reported in the text. Whereas fa, ji, and are obtained by means of 
scaling laws, using the most accurate theoretical estimates for standard critical exponents [54] 



with our Pade final results of the RG dimension y» 





n = 


n = 1 


n = 2 


n = 3 


n = 4 


n = 5 


n = 16 


V2 (PBL) 




1.7348(2) 


1.7645(3) 


1.7897(3) 


1.8112(6) 


1.830(1) 




3/2 (Pade) 




1.733(6) 


1.763(4) 


1.789(3) 


1.811(2) 


1.830(2) 


1.927(4) 


4>2 


1 


1.092(4) 


1.184(3) 


1.272(2) 


1.356(4) 


1.398(4) 


1.755(4) 


02 




0.798(4) 


0.831(3) 


0.861(2) 


0.891(3) 


0.894(3) 


0.978(4) 


72 




0.294(8) 


0.353(5) 


0.411(4) 


0.466(3) 


0.504(3) 


0.778(7) 


2/3 (PBL) 


0.7258(2) 


0.8153(3) 


0.8920(2) 


0.956(2) 


1.011(2) 


1.059(4) 




2/3 (Pade) 


0.725(30) 


0.814(10) 


0.891(9) 


0.957(4) 


1.014(4) 


1.063(7) 


1.31(1) 


03 


0.426(18) 


0.513(6) 


0.598(6) 


0.681(3) 


0.759(4) 


0.812(6) 


1.19(1) 


As 


1.336(18) 


1.377(6) 


1.416(6) 


1.453(3) 


1.488(5) 


1.480(7) 


1.54(1) 


73 


-0.91(4) 


-0.865(13) 


-0.818(12) 


-0.772(6) 


-0.728(6) 


-0.668(11) 


-0.35(2) 


j/4 safe 


-0.380(18) 


-0.23(1) 


-0.098(6) 


0.012(6) 


0.104(8) 


0.188(8) 


0.62(1) 


j/4 best 


-0.393(5) 


-0.236(3) 


-0.103(1) 


0.0094(30) 


0.107(5) 







, 328835. +90677.3 n+10871.8 n 2 +1595.82 n 3 +446.805 n 4 +28.9975 n 5 +0.574653 n 6 ^4 

(n+8) 7 T 



+ 
+ 



-3.06135 10 7 -1. 35901 10 7 n-2. 55849 10 6 n 2 -179042. n 3 +24766.4 n 4 +8290.96 n 5 +672.157 n 6 +22.2622 n 7 +0.318104 n 8 5 

(n+8) 9 ' 

3.18214 10 9 +1. 83869 10 9 n+4.52163 10 8 n 2 +5.94633 10 7 n 3 +5. 79335 10 6 n 4 +913345. n 5 +161861. ra 6 +14567.2 n 7 

(n+8) 11 

+669.512 n 8 +17.0980 n 9 +0.199456 n 10 " 
(n+8) 11 



(10) 



From a general analysis [1,8], it is known that 2/4 is positive for n > N c and negative in the 
opposite case, with iV c a bit smaller than 3 [1]. The series (10) is alternating in signs for 
n < 5. However, in this case, the coefficients are not so small for a simple Pade summation 
to be effective. Thus, we apply a PBL resummation. Fig. 1 sketches the non-defective 
PBL approximants for n = 2, 3. It is evident that several approximants are very close 
to each other, whereas the [1/M] ones are well separated. Usually, when facing with a 
similar situation, the PBL approximants that are far from the mean-value are discarded in 
the averaging procedure. Since in this case we may not take advantage of the Pade data, 
in order to be sure to not underestimate the error, we report in Table IV two different 
estimates: one is the average over all PBL (which we will term "safe"), the other ("best") 
is obtained discarding the [1/M] approximants. For n — 5, 16 we report only the "safe" 
estimate, because we do not have anymore evidence for a trend in the [1/M] approximants. 



IV. CONCLUSIONS 

In this paper we determined the critical exponents associated with harmonic operators 
of degree 2, 3, and 4 in O(n) models by means of pseudo-e expansion. All our results for 
fa, 7i, and $ are reported for n = 0, 1, 2, 3, 4, 5, and 16 in Table IV. 

In order to make a comparison with the values in the literature we report in Table V 
all the most accurate theoretical estimates for yi, as obtained by means of scaling relations 
(4) using the most precise determinations of standard critical exponents [54]. For all values 
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TABLE V. Theoretical estimates of the RG dimension j/j as obtained by various approaches for 
several n: five-loop e expansion (e exp), six-loop fixed-dimension expansion (FD), high-temperature 
expansion (HT exp), Monte Carlo simulations (MC), and 1/n expansion at order 0(l/n 2 ) or 
0(l/n). The results are obtained by using scaling relations, using the most precise theoretical 



estimates for standard critical exponents [54]. 
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l.olo(b ) 


1.8oz(8J 




IVLO [4LZJ 
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MC [44] 












1.815(39) 




HT exp. [9] 






1.750(22) 


1.758(21) 








0(l/n) [57] 








1.640 


1.730 


1.784 


1.932 


0(l/n 2 ) [46,58] 








1.78(14) 


1.83(6) 


1.88(2) 


1.95(3) 


2/3 


n = 


n = 1 


n = 2 


n = 3 


71 = 4 


n = 5 


n = 16 


6-loop (pscudo-e) 


0.725(30) 


0.814(10) 


0.891(9) 


0.957(4) 


1.014(4) 


1.063(7) 


1.31(1) 


6-loop (FD) [34] 


0.758(19) 




0.895(15) 


0.953(23) 


1.015(31) 


1.065(19) 


1.310(13) 


5-loop (e-cxp) [34] 


0.739(9) 




0.892(22) 


0.958(42) 


1.020(45) 


1.064(25) 


1.28(10) 


0(l/n) [57] 










0.791 


0.933 


1.323 


2/4 


n = 


n = 1 


n = 2 


n = 3 


n = 4 


n = 5 


n = 16 


6-loop (pscudo-e) best 


-0.393(5) 


-0.236(3) 


-0.103(1) 


0.0094(30) 


0.107(5) 






6-loop (pscudo-e) safe 


-0.380(18) 


-0.23(1) 


-0.098(6) 


0.012(6) 


0.104(8) 


0.188(8) 


0.62(1) 


6-loop (FD) [56,55] 






-0.103(8) 


0.013(6) 


0.111(4) 


0.189(10) 




5-loop (e-exp) [56,8] 






-0.114(4) 


0.003(4) 


0.105(6) 


0.198(11) 




MC Ref. [43] 






-0.17(2) 


-0.0007(29) 


0.130(24) 






0(l/n) [57] 












-0.08 


0.662 



of n our estimates are in perfect agreement with all known results and in the majority of 
the cases they are the most precise ones. We stress that such high accurateness should 
not be due to underestimation of the uncertainty, since we check our final error bars with 
other resummation techniques, such as Pade-Borel-Leroy and conformal mapping. However 
the large n results need to be discussed, since our estimates for n ^ 6 can be affected by 
systematic errors, because the perturbative series we summed have not alternating signs. 
Anyway, for n — 16 all the theoretical estimates are in good agreement, signaling that the 
evaluation of our uncertainty is probably good even in this case. 

Let us finally compare our values with some experiments. We mention the result 2 =2 — 
1.17(2) for the (2—^1 + 1) bicritical point in GdA10 3 [59], and = 1.279(31) in the 
(3 — > 2 + 1) study of MnF 2 [60]. Other experimental measures of 2 can be found in Ref. 
[61]. The experimental results obtained for a nematic-smectic-A transition reported in Ref. 
[22] are /3 2 =2 = 0.76(4) and 7 2 =2 = 0.41(9). For the third harmonic exponent we quote 
/3™= 2 ~ 1.66 in liquid crystals [23], (3%= 2 = 1.50(4) (Ref. [32]) (3%= 2 = 1.80(5) (Ref. [30]) in 
Rb 2 ZnCLi. All these values compare well (within their own uncertainties) with our results. 
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